Interedge coherent line junctions in Quantum Hall systems 
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In this paper we address some of the properties of quantum Hall line junctions (QHLJ) that occur near barriers 
separating electron gases on quantum Hall plateaus. In narrow barriers where electron tunneling can occur, the 
low energy physics of the QHLJ is described by the quantum sine-Gordon model. We propose procedures to 
study a sort of properties of these systems in relation with recent experimental studies of Kang et al, Nature 
403, 59 (2000). We propose experimental ways of measuring the values of the dimensionless coupling constant 
characteristic for the Sine-Gordon model as well as its Fermi velocity. When the bulk filling of the 2DEG 
subsystems is v ~ 2, these systems can be good candidates for observation of spin-incoherent LL behavior. 

PACS numbers: 73.43 Jn, 73.43.Cd, 73.43.Lp, 71.10. Pm 
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I. INTRODUCTION 

The edge states of quantum Hall systems offer a rich ground 
for testing the theoretical predictions on the properties in one 
spatial dimension (ID) of strongly interacting fermionic sys- 
tems. Unlike in quantum wires on which impurities can lead 
to localization effects hindering the observation of the generic 
properties of ID fermions, the QH edge states are considered 
to be very clean realizations of the Luttinger liquid descrip- 
tion. The QHLJ that occur near disorder free barriers sepa- 
rating electron systems on quantum Hall (QH) plateaus, have 
been shown to be an important tool on these studies. 

Properties of QHLJ, as the differential tunneling conduc- 
tance between counter-propagating edge states has been stud- 
ied extensively in experiment by Kang and coworkers 1,2 , 
by Grayson et al.r^, and by Roddaro et air- as well as 
in theor y 6 ' 7 ' 8 ! 9 . In the experiments of Refs. OB the two- 
dimensional electron gases (2DEG) are separated by a nar- 
row barrier, constructed with the atomically precise cleaved- 
edge overgrowth (CEO) technique. Such narrow barriers en- 
able strong correlations between electrons on opposite sides 
of the junction. At finite barrier heights the intersections of 
the Landau levels in the barrier region turn into gaps. When 
v falls within the gap, transport along the barrier is forbid- 
den and the QH effect extends throughout the sample. In 
Ref. Ill the tunneling conductance across the barrier displays 
zero bias peaks of amplitude G ~ 0.le 2 /h over a range of 
filling fractions 8v w 0.3 centered at v* = 1.35 [as well as 
at points between other higher integer fillings]. The size of 
the gap (A ~ 8v huj c ) however has been a point for which the 
theoretical 7,8 and experimental 1 results have been in partial 
disagreement. 

In Ref. [ZJ the formation of a phase-coherent state of 
fermions between opposite sides of the barrier has been pro- 
posed for the explanation of the I-V characteristics of this 
system. According to this model, in the presence of fermionic 
tunneling between edges on opposite sides of the barrier, the 
low energy physics of the QHLJ is expected to be described 
by the quantum sine-Gordon (SG) model in terms of variable 
ip that represents the phase of the interedge phase-coherence 
order parameter. The properties of the SG model are depen- 
dent on the value of the coupling constant [3 and the value of 
the Fermi velocity. The accuracy of estimation of the physical 



mass in this quantum field theory however suffers from two 
limitations. First, the SG parameters (3 and vp appearing in 
the bosonized version that describes the line junction at low 
energies are hard to be estimated with accuracy from micro- 
scopic models. Second, the physical mass depends also on the 
choice of the ultraviolet cutoff regulating the short distance 
behavior of the quantum field theory. Progress in the latter 
direction has been achieved by using the method of asymp- 
totic matching 8 , in which the long distance of the quantum 
field theory and the short distance of the microscopic model 
are matched in the region of common validity. The theoreti- 
cal calculations for the magnitude of the gap within this model 
give an estimate of A - 1.3K (Ref.0) and A ~ 1.5K (Ref.H). 
In the latter calculation, the increased value of magnitude of 
the gap was attributed to the renormalization of the tunneling 
amplitude by interactions as well as to the better estimation 
of the short-distance cutoff of the model. The experimental 
value of the energy gap above the ground state measured in 
the experiments of Ref. [l| is one order of magnitude higher 
A ~ 15K. To explain the features of the experiment, in Ref. 9 
the model with single strong tunneling center along the barrier 
was proposed. Wide barriers due to differences in the Fermi 
wave-vectors of the states on opposite sides of the barrier, do 
not allow for tunneling to take place at any point except at 
imperfection ones created unavoidably during the fabrication 
process. Narrow barriers on the other hand can allow for such 
tunneling to take place at any point along the barrier. Our 
work here is intended for the case when tunneling can take 
place at an infinite number of points along the barrier. The 
results and predictions of this work if examined in the experi- 
ments can however shed light on the model that best describes 
the experiments. 

Assuming therefore the model where tunneling can take 
place at an infinite number of points along the barrier here 
we show that one can indirectly extract information on the 
value of these quantities (J3 and vf) by using the integrability 
property of the SG model. The spectrum of the quantum sine- 
Gordon model consists of massive topological particles, soli- 
tons, antisolitons, and when forward interactions are strong 
enough also of their bound states. In the presence of a chemi- 
cal potential however that couples with these charges the spec- 
tral gap can be suppressed. This phenomenon is similar to the 
one that is discussed previously in the context of quantum fer- 
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romagnetism and phase transitions in bilayer QH systems in 
the presence of a magnetic field coplanar with the 2DEG lay- 
ers and also interlayer tunneling experiments in bilayer QH 
systems in the presence of a magnetic field coplanar with the 
2DEG layers (see Refs. llOlll II) . When this chemical potential 
exceeds the mass of the soliton (A /2), a finite density of soli- 
tons appears in the ground state, distributed on a Fermi sea 
according to their statistics, embedded in their interactions (or 
their scattering matrix). The low-energy physics of this sys- 
tem then will be of particle-hole type formed around these 
Fermi points. The properties of this metallic state, namely 
the value of the Luttinger liquid (LL) parameter K and the 
Fermi velocity can be accessed with the thermodynamic Bethe 
ansatz. It turns out that at values of chemical potential not very 
far from A, the LL parameters K and vf approach the values 
of the SG coupling constant (3 and vp prior to gap opening. 
Experimentally there are two quantities that offer the measure- 
ment of two combinations, the product and the ratio of the LL 
parameter K and the Fermi velocity, namely the Drude weight 
and charge susceptibility, respectively. 

On the other hand, in the insulating phase of the model we 
propose the measurement of the optical conductivity of the 
system. The optical conductivity measurements should mea- 
sure a delta function-like presence for each of the kinds of 
particles of the spectrum odd under charge conjugation and 
also a contribution on the form of a continuous band for their 
particle-particle combinations (with specific selection rules). 
The proposed measurement of the optical conductivity there- 
fore would be an important confirmation of the existence of 
the interedge coherent state. It would serve also the purpose 
of bounding the value of the coupling constant. 

The SG model has applications in other contexts too where 
therefore it would be helpful to know when the finite tem- 
perature LL description starts to become unreliable. This is 
something that we also discuss in this paper. 

At filling fractions v ~ 2 the spin degrees of freedom play a 
role in the dynamics of the system. In this case terms like spin 
exchange interactions as well as the coupling of the spin cur- 
rent to the magnetic field lead to a rich physics. As discussed 
in Ref. 9, the spin sector can be in a gapped state. This gap can 
be suppressed in the presence of the magnetic field. At values 
of this field only slightly in excess of the gap the energy scales 
in the spin and the charge sectors separate. At a small window 
of energies above the spin gap at finite temperatures the spin 
incoherent LL (SILL) state 12 , can be realized. Experimental 
implications on charge and spin transport are discussed. 



II. THE MODEL 

In this section we introduce the model that intends to de- 
scribe the experimental systems used to study properties of 
QHLJ. The geometry of the model is represented schemati- 
cally in Fig. Q] A QH system at a plateau at integer filling 
fraction v = 1 is separated down the middle by a narrow bar- 
rier as illustrated in Fig. Q] The barrier creates on its opposite 
sides edge states that propagate in opposite directions. The 
potential that CEO barriers present to the electron gas in the 



experimental systems are sufficiently abrupt that edge recon- 
struction phenomena and complications in their description 
do not arise. Therefore we describe the edges by single chi- 
ral moving channels. At the integer filling fraction v = 1 the 
spin degree of freedom is frozen and does not contribute in the 
dynamics of the system at low energies. 

The QH samples used in Ref. 1 are long in the direction 
of the barrier. More specifically, the lateral dimensions of the 
QH samples used are ~ lOO^m long in the direction of the 
line junction (noted L x below) and 13— 14/um wide from each 
side of the junction. The 2D electron systems are separated 
from each-other by a 8.8nm thick barrier. 

The electron states on each side of the barrier can be 
described by single-particle wave functions extended in the 
direction of the barrier - 0fe( r ) = e ikxX Xk(r)/ 'y/Lx, chosen 
here to be the x direction, and labeled by an one-dimensional 
wave vector k x that is proportional to the guiding center 
along which the wave function's y-coordinate is localized, 
Y = k x l 2 . In mean field Hartee-Fock theory the dispersion 
relation of electrons can be found by solving the equation 

{-^l + ^(y - l 2 kf + V H {y) + V B {y)} Xk (y) 

Jdy'V E (y,y')xk(y') = e k \k(y), where V E (y,y') is the 
exchange interaction experienced by the electrons on the 
same side of the barrier (described by the Slater determi- 
nant). where V E (y, y') = (2/L x )J2' k xi,(y')Xk>(y)K(k- 
k)'(y — y')) is the exchange interaction experienced by 
the electrons on the same side of the barrier and V H (y) = 

(-l/£*)E fe Jdy'|xH2/)| 2 M((y - V'f + ™ 2 )/L x ] is 
the direct interaction. Vb{v) is the potential presented to 
the electrons by the barrier, whereas w is a short distance 
cutoff of the order of the magnetic length. The dispersion 
relations e{k x ) and the corresponding gaps that appear in the 
intermixing region, are shown in Fig.Q] right. 

The classical quadratic fluctuations of the coupled chiral 
edges on the sides of the barrier are described by a Hamil- 
tonian density of the form H ~ u a {x)K a p(x — x') up(x'), 
where summation for repeated indices is assumed. u± are the 
transverse displacements of the right and left edge, respec- 
tively, from the ground state location. In principle these inter- 
action kernels can be calculated by using of the Hartree-Fock 
mean field theory. 

Adopting field theory terminology the quadratic Hamil- 
tonian can be written in terms of chiral currents J a (x) = 
u a (x)/(2irl 2 ), where a = ±1, for the current on the right 
and the left of the junction, respectively 

H= fdxdx' V J a (x)K a0 J p (x') . (1) 

J a,/J=± 

The Hamiltonian can be rewritten in terms of sums and differ- 
ences of currents of the edges on the left and on the right of 
the barrier, with the corresponding kernels also given by sums 
and differences of the kernels of same and different sides, 
K 1 (x — x') = l 2 [K RR (x ~ x') + "/K RL (x — x')], where 
7 = ±. Since the barrier's thickness is of the order of the mag- 
netic length, the Coulomb interaction between fermions on the 
same side of the barrier K RR and Kll, where K RR = Kll, 
are comparable to the one between fermions on opposite sides 
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FIG. 1: Left: Schematic illustration of a QHLJ formed in QH 
systems by means of cleaved edge overgrowth barriers, similar to 
the ones used in experimental studies of Kang and coworkers 1 . 
The barrier is considered to be narrow enough to allow for in- 
teredge tunneling to take place at an infinite number of points. 
Right: the electron dispersion curves e(k x ) in the presence of 
interactions are shown. In mean field Hartee-Fock theory the 
dispersion relation of electrons can be found by solving the equa- 

tion {-*Ldl + ^{y-l*ky + V H {y) + V B (3j)}xk{y) - 

fdy'V E (y,y') X k(y') = e kX k(y), where V E (y,y') = 
(2/L x )J2' k xl>(y')Xk'(y)K-(k - k')(y - y')) is the exchange in- 
teraction experienced by the electrons on the same side of the barrier, 

v H ( y ) = / dy' \xAyTH{{y-y') 2 +™ 2 )lLl], 

is the direct interaction, and Vb (y) is the potential presented to the 
electrons by the barrier, w is a short distance cutoff of the order of 
the magnetic length. At finite barrier height, energy gaps appear in 
the barrier region. 



of the barrier Krl. Therefore one of the kernels K— (Kq 
below) should be expected to be much softer than the other. 
The Hamiltonian (HJ is quantized by requiring for the chiral 
currents to fulfill the commutation relations 



[J a (x),Mx')} 



a—6 a pd x 5{x ~ x') 
Ztt 



(2) 



The chiral currents J(x) =: W(x)R(x) : and J(x) —: 
L'(x)L(x) :, where W and are the fermion creation op- 
erators on the right and left moving edge, respectively, can be 
expressed in terms of bosonic fields 



J(x) 



-.d x (j) R (x) , J(x) 



-.dMx) . (3) 



Here we used the convention R(x) ~ e -iv / 47r0 f j(a;) f or ^ 
right movers, and L{x) ~ e %V^<t>L{x) f or t jj e i e f t movers 
(up to factors of l/y/2n) where 4>r(x) and <Pl(x) are the 
chiral components of the non-chiral bosonic field <!>(x) = 
4>r( x ) + <Pl(x) that is formed along the barrier. Intro- 
ducing the dual bosonic field Q(x) — 4>r(x) — 4>l(x), 
which satisfies the commutation relation [O(x), d x Q{x')\ — 
—iirS(x — x'), the Hamiltonian density can be written as a 
sum H ~ d x $(x)K<s,(x - x')d x ><S>(x') + d x S(x)K e {x - 
x')d x 'Q(x'). The collective modes of the system are co(q x ) = 

±Trl 2 [K< s >(q x )Ke(qx)} 1/2 \qxl where K$(q x ) and Ke(q x ) 
are Fourier transforms of the interaction kernels. The Fermi 
velocity in the case of long range interactions is momentum 
dependent. 



To make some theoretical progress we assume here never- 
theless that the interactions are short ranged. Such an assump- 
tion should not lead to significant changes on the physics of 
the system. Ti in this case takes the usual LL form with mo- 
mentum independent Fermi velocity hvp = (K^Kq) 1 / 2 and 
a LL constant K = (K / 'K®) 1 ' 2 . 

The Hamiltonian describing the system along the barrier in 
the the presence of tunneling H t = toR'L + h.c. is given by 

H = J dx [(d x <S>) 2 + {d x e) 2 ]+2ficos(V^K<i>(x)) , 

(4) 

where /i = f(to,K). The cos operator extends along the 
whole length of the barrier. Only at the noninteracting point of 
K = 1, is /i = io/27r. The exact functional dependence \i = 
f(to, K), as shown by (Ref.|8|), is dependent on microscopic 
details of the underlying system. We estimate the value of the 
LL parameter if to be given by AT 2 ~ (d 2 /I 2 )/ ln(L x /l), and 
differences in the value of // from the free point value of to /27r 
are expected. 

Eq. (0| represents the well-known sine-Gordon model. In 
this model the relevance of the electron tunneling processes 
between edges on opposite sides of the barrier depend on the 
value of the parameter K. This parameter on the other hand 
depends on the value of the filling fraction of the QH liquids 
on the sides of the barrier, the interaction strength along the 
QH line junction, and the interaction strength between charges 
of counterpropagating edges on each QH subsystem on each 
side of the barrier. From previous calculations 13 the interac- 
tions across the Hall bar and the interactions along the line 
junction have the opposite effect on the relevance of the tun- 
neling processes. In general there can be geometries in which 
these interactions can be tuned to cancel each others effect. 
It is assumed however that the Coulomb interactions between 
charges located on the edges on the top and bottom of the Hall 
bar in systems examined by Kang et cd x , due to their spatial 
separation, are small thus can be neglected in the following. 

The sine-Gordon model obtained above is a well-known 
model in quantum field theory and many of its properties are 
known. This model has a discrete symmetry $ — > $ + 27rn//3, 
where 1 — 4ttK. At f3 2 /8ir < 1 in the ground state the 
symmetry is spontaneously broken and the spectrum becomes 
massive. In the following we absorb the factor 1/8tt in the 
definition for f3. In the region 1/2 < 1 < 1, known as the 
"repulsive" regime, the spectrum is composed of topologically 
charged particles called solitons and antisolitons. In the region 
< f3 2 < 1/2, known as the "attractive" regime, reached in 
the presence of forward scattering interaction, the spectrum is 
enriched by the their bound states, called breathers, the num- 
ber of which increases with decreasing (3. The mass spectrum 
of the model is 



M n = 2M s sin ri- 



ft 2 



2 1 - p 



, n= 1,. 



1-/3 2 

p 2 



(5) 



where n is integer and [•] represents the integer part of a real 
number. With M n we denote the mass of the n-th breather, 
whereas with M s we denote the mass of the soliton. The soli- 
ton and antisoliton have equal mass. 
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The scaling dimensions of the cos operator are d = K/2 
and at K < 2 the QHLJ turns into an insulator along its 
length when the chemical potential for topological particles 
falls withing the gap. In this case the QH effect is established 
along the whole sample and current flows along its perimeter. 
Therefore at zero temperature one would expect the height of 
the zero-bias conductance peak to be e 2 /h. At finite temper- 
atures however there will be some particles in the spectrum 
even at chemical potential below the mass gap. As a result 
even when the chemical potential falls within the gap, current 
will partly be carried along the junction. When the chemi- 
cal potential falls outside the gap the QHLJ becomes an ID 
conductor and at low voltages and low temperatures the con- 
ductance across the barrier vanishes. 

The mass of the physical particles of the SG model, in the 
context of the QHLJ systems with interedge tunneling present, 
constitutes half of the gap that is measured in the experiments 
of Ref . 1 . It depends on the microscopic parameters of the sys- 
tem, the exact value of the Fermi velocity, coupling constant 
and tunneling amplitude as follows^ 




where A is a short distance cutoff of the field theory. These 
quantities however are very hard to estimate with accuracy 
in simplified microscopic models. This is the reason that we 
try to present here indirect ways of estimating these quantities 
experimentally. 



III. EXPERIMENTAL POSSIBILITIES FOR ESTIMATION 

OF (i 2 ,v F 

As it was discussed in Refs. [Hll even though the SG 
model has a massive spectrum, massless modes can be in- 
duced in the presence of external fields that couple to the 
charges of the topological excitations. We show below how 
the parameters (3 and vf can be extracted using the thermody- 
namic Bethe ansatz method. In QHLJ that we are considering 
here, we discuss two situations where the gap suppression oc- 
curs. 

I. The first one is the case when the filling fractions on the 
2DEG subsystems is v = 1. In this case such a field can be 
some linear external charge Q(x) that couples to the charge of 
the ID non-chiral LL that extends along the separating barrier, 
and can be represented as 

« a = /d,d.W(*-*W) • (7) 

V(x — x') is an interaction kernel. Just as in the case of Eq. © 
we assume here that the interactions are short ranged due to 
the presence of screening gates or other metals in the vicin- 
ity of the barrier. The total charge density as usual in the 
bosonization approach is represented by 

p(x) = po + -^=d x $(x) , (8) 



where the constant charge density po = hp /it is usually neu- 
tralized by the positive background charge. 

II. The other possibility that we will consider below arises 
at bulk fillings v ~ 2 of the 2DEG subsystems. In this case 
the electron spin is not frozen as in the previous case. The spin 
degree of freedom plays an important role on the dynamics of 
the edges. The sine-Gordon model can arise here either in the 
presence of antiferromagnetic interactions or in the presence 
of ferromagnetic type interactions with magnetic anisotropy 
(Ref. |9J), as we elaborate below. The external field that can 
drive the transition to the gapless phase can be the magnetic 
field (that couples with the spin current) 

H Zccm = -fi B gB(J z + ,!-) = ^^-d x ^ s (x) , (9) 

v 7r 

where J z + J z = — d x $ s /y/ir was used. In Eq. (O g e is the 
gyromagnetic ratio of the electron and /is is the Bohr magne- 
ton. 

For the field that couples with the gradient of the field d x &, 
whether we will be talking for the first or second case, we will 
use the notation H, where we will mean either H ~ QV or 

H ~ PB9eB. 

The presence of the H field in the Hamiltonian, has im- 
portant consequences in the physics of the system. In the 
absence of H, the ground state of the system is reached at 
the trivial, the static and constant in space field configura- 
tions Q(x) — 2irn/(3, being the discrete degenerate minima 
of the potential energy. The presence, on the other hand, of 
Hd x &(x) and the harmonic term in the Hamiltonian as op- 
timal require the new field configuration be <£>(a;) ~ Hx. 
These competing tendencies are resolved at a critical value 
H c of H, above which the system is characterized by a fi- 
nite value (d x Q(x)) ^ (or finite density of solitons). In the 
new phase the field $ generally follows the Hx line, but such 
as to minimize also the potential energy cos (3Q(x), leading 
to finite regions of peaks in the charge density. This is the 
commensurate-incommensurate transition, first discussed in 
Refd 

On the quantum level, H serves as chemical potential for 
the topological particles. In its presence the solitons lower 
their energy, and above a critical value of the chemical poten- 
tial, they will tend to proliferate in the ground state. These 
particles however interact with each other, leading to a Fermi- 
sea-like distribution in the ground state. Gapless modes would 
consist of excitations above this Fermi sea. Our aim here is to 
find the value of the parameters of the SG model, (3 and vp, 
discussed in the previous sections. With the thermodynamic 
Bethe ansatz, we can calculate the values of the LL parameter 
K of this metallic phase, as well as the corresponding Fermi 
velocity vp, as functions of both, chemical potential as well as 
the SG parameters (3 and vp. K and vp at large values of the 
chemical potential relative to the mass gap approach respec- 
tively values of (3 2 and vf of the insulating phase. However, 
we show that already at H ~ 1.5(A/2) these quantities are 
practically to within 5% range of the values of (3 2 and vp of 
the SG model. 

From the practical point of view, in the experimental stud- 
ies, one way of getting access to the LL parameters in the 
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H=0.1M, 




FIG. 2: On the left, the energy spectrum of the SG model in the 
presence of a small external field H, here H/M s — 0.1, to allow 
for a small split in the spectrum, is shown. In the absence of an 
external field and at T = the soliton and antisoliton energies coin- 
cide, whereas the breathers energy is slightly larger [see Eq. ((5)]. For 
this representation the temperature was taken to be T/M a = 0.01. 
On the right, the spectrum in the presence of an external field that 
exceeds the soliton mass, is shown. For H > M 3 the soliton gap 
closes and in the ground state a soliton condensate forms. Due to the 
opposite sign of their topological charge the energy of the antisoli- 
tons increases. The system is at T — or very small T, compared 
to M B , a solitonic Luttinger liquid. The graphs are obtained from the 
numerical solutions of the defining integral equations. 



incommensurate (I) metallic phase is by measuring the Drude 
weight, which is given by the product vpK, Another way 
to access a different combination of v F and K is by mea- 
suring the charge susceptibility, which is given by their ratio 

X~K/v F . 

In the following first we discuss the applicability of the LL 
picture at finite temperatures. 



A. Luttinger Liquid picture of sine-Gordon model in the 
presence of external fields: Finite temperatures 

Before discussing the question of the Drude weight and 
charge susceptibility in the metallic phase of the SG model we 
first discuss the applicability of the LL picture at finite temper- 
atures. The case that we are discussing, the LL is obtained by 
linearizing the spectrum around the two Fermi points of the 
soliton sea. This approximation works at its best at T = 0, 
and its accuracy decreases with increasing temperature. The 
experiments from the other side are never performed at zero 
temperatures. It is useful therefore to examine its applicability 
at finite T. 

With increasing temperatures the particles start filling states 
away from the Fermi surface and at relatively higher temper- 
atures the breathers and antisoliton branches start filling too. 
As a result the low energy excitations are not of a particle-hole 
type close to the Fermi surface any more. Here, the applica- 
bility of the LL picture is examined by studying the depen- 
dence of the susceptibility \ on T. Its divergence at a certain 
critical value of the external field notifies a transition from a 



gapped system to a gapless one with sharp Fermi surface. At 
finite temperatures there will be some particles in the spec- 
trum even at external fields below the mass gap, as a result 
the divergence turns to a finite height peak which decreases 
continuously with temperature. 

To obtain the dependence of the susceptibility on the exter- 
nal field H and temperature T exactly, we use the thermody- 
namic Bethe Ansatz (TBA), a summary of which is contained 
in the Appendix. 

Here we examine the SG model only at the reflectionless 
points f3 2 = 1/n, where the particle scattering is diagonal 
(except the asymptotic case H/T ^> 1, which is applicable 
for any (3). In the following the calculations are performed in 
rapidity space 9 which parametrizes the energy and momenta 
as follows 



E a {6) = m a cosh< 



VF 



preserving thus the relativistic dispersion relation 

E 2 a (P a )=v F P 2 + ml . 



(10) 



(11) 



m a are particles' masses and vp their velocity. Before pro- 
ceeding any further a few words are in order on the equations 
for the spectral functions obtained with the Bethe ansatz tech- 
nique that we will use in the following. The Bethe ansatz 
equations are obtained from the set of periodic boundary con- 
ditions set on the phases accumulated by particles with dif- 
ferent momentum as they travel around the circumference of 
a circle of length L in the presence of the other particles. In 
these equations the interactions are reflected in the phase col- 
lected as two particles interchange positions. First of all, it 
is a general fact that no two particles in integrable theories 
can occupy the same state (for then the Bethe wave function 
that describes the overall system vanishes). Therefore in these 
models, in the ground state, the particles are distributed on a 
Fermi sea. The energy spectrum is determined by a complete 
system of coupled integral equations 



rn. 



cosh6> - H a q a = e a 



1 



-t b /T 



(12) 

where (*) is a short notation for the convolution of K a b(0—6') 
with In [1 + exp{— e&(0')/T}]. In Eq. (O the subscripts a, 
b stand for s, s, for soliton, antisoliton, and n for the rt-th 
breather, respectively. H a is a chemical potential for the parti- 
cles of the spectrum. If interactions are neglected the spectral 
function agrees with the bare one [left-hand side of Eq. (fT2bl. 
Their topological charges are q s = —q s = 1 and q n = 0. 
Due to the sign difference of the charges of solitons and an- 
tisolitons in the presence of a chemical potential their spectra 
become different. Particularly at H > M s , the energy spec- 
trum will contain an interval of states of negative energy. At 
T = and H > M s , due to the fact that only one of the spec- 
tra takes negative values between the points ±B (the Fermi 
points), the equations (fT2l simplify. The term involving K a {, 
on the right-hand side of Eq. ( fT2T > reduces to [K ss * e s ](6) in 
the interval (—B,B). The Fermi points are defined self con- 
sistently from e s (±B) = 0. The Kernels K a b(9) are obtained 
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from the scattering matrices by taking 

K ab {9) = -^-^\nS ab {9) . (13) 

Z7T U(7 

The scattering matrices encode all the information about the 
particle interactions. For the sine-Gordon model they are 
known (Refs. 18.19 20 21), and for completeness we give 
them in the following 



2tt 



[S n 



K n ,± = - 



cosh (£§cj) cosh[(l - <)§ w] sinh (Z£f c 
cosh (f sinh (£f w) 

dwe l " e cosh (£f w) sinh (n£§ w) 



2tt 



cosh (f u; 



sinh (ff w) 



(14) 



A" 



dme^ 9 sinh [(1- Qf a;] 
2tt 2 cosh (fa;) sinh 



1; n > I 



ft 2 



1-I3 2 ' 

In Eqs. d3"5ll36i > below, the minus sign is absorbed in the Ker- 
nel. 

The Kernels can be calculated easily for the cases when 
(3 2 = 1/ n, for n integer. In the Appendix we give the Kernels 
for the n = 3 and n = 5 when the spectrum contains just one 
and three breathers, respectively. 

The integral equations for the energies of the particles 
Eq. ( fT2l can be trivially solved at the free fermion Luther- 
Emery point (3 2 = 1/2. At all other values of (3 2 they can be 
solved numerically. The number of coupled integral equations 
increases increases with decreasing value of the coupling con- 
stant f3 2 . 

In Fig. [3] (left) we give plots of the energy spectra at (3 2 = 
1/3 and k B T = eV = 10M S , where V = 2H. In Fig. [3] 
(right) also the particles' distribution functions are shown. 
The fact that no two particles with same momenta occupy the 
same state is an well-known feature of the integrable models, 
resembling this way the fermionic statistics. In fact they have 
a mixed or anyonic statistics since their scattering matrix de- 
pends on the difference of their rapidities. These particles are 
neither fermions nor bosons. The phase accumulated by in- 
terchanging these particles will depend on the difference of 
their rapidities. Notice, for instance, that interchanging parti- 
cles with 61 — 62 — > +00, results in Sq(9i — 62) = 1, thus 
resembling bosons. 

From BA, for the susceptibility we have 



X = 



d 2 f(p 2 ,H,T) 
dH 2 



(15) 



where / is the free energy of the system in the presence of a 
chemical potential. In TBA it is given by 



+ OO 



d0 

27 



cosh 9 In 



1 



(16) 



T=10, V=10, M_s=1 

\ 
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FIG. 3: On the left panel the energy spectra of solitons, breathers 
and antisolitons, respectively, in rapidity space are presented. On the 
right panel the respective distribution functions are presented. The 
particles behave as fermions for 9 — » 0, no two particles of the same 
kind with the same rapidity can exist. These particles are neither 
fermions nor bosons. The phase accumulated by interchanging two 
particles depends on the difference of their rapidities. At zero tem- 
perature and in the presence of field V > 2M B the antisoloton and 
breather branches remain unpopulated whereas the solitons reside 
in the interval (— B, B) with a step function distribution in rapidity 
space. B is defined by the equation e 3 (±B) — and Eq. d 1 2b for 
T = 0. At T 2> M s , which is the case presented here, the antisoliton 
and breather particles are also present. 



where e a (9) is the energy of the particles in the presence of 
interactions. 



symptotic limit H/M s 3> 1 of charge susceptibility at T 7^ 



1. A 



In the limit of large values of chemical potential above the 
mass gap, the functional dependence of the free energy (and 
therefore of the susceptibility) on (3 2 , and H, can be obtained 
analytically. 

In the H/M s ^> 1 limit and temperatures of the order of 
the soliton mass, T ~ M s , one can assume that the spec- 
trum contains only solitons. Breathers and antisolitons are 
also present, but as can be inferred from Figs. [2] [3] their 
density is small due to high energies required to fill these 
branches. The formula for the free energy Eq. dT2b in this 
limit simplifies to a single integral, accounting simply for the 
energy of the solitons. Moreover, due to the large density 
of states available to solitons at higher rapidities (reached at 
high H), the main contribution in the free energy comes from 
the region cosh 8 ~ H/M s — 1, Here we can approximate 
cosh 9 ~ e e /2, and rewrite the equation for the soliton energy 



e s (9)+TK ss *ln 



1 + e 



(0) 



M s e & 



H. (17) 



In this limit it is assumed that all solitons have v ! F 



v F , 



and the TBA equations and e s (6), resemble the massless case 
of Fendley et al*^. This e s (0) differs from the original one 
at the low energies. It turns out that e s (— 00), will enter the 
expression for the functional dependence of / on H, (3 2 , and 
T. 

Using Eq. (1171 , we can split the free energy into two parts. 
After technical steps 30 , we find that at H/T ^> 1, the first one 
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/i converges to the H independent constant fx = —irT 2 /12. 
In fact the leading order terms in fi are 



h 



ttT 2 
12 



T 2 

£_ e s (-oo)/T. 

2tt 



re »(- 00 ) c c(-ooVr , 

47T 



(18) 

The second part, with H contained in leading order, turns out 
to converge to fa — —0 2 H 2 /tt. This is the part that we are 
interested here. More precisely f% is given by 



h 



HT 
~2V 



In 



1 



,-e,(-oo)/T 



(19) 



where e s (— oo) can be found from Eq. dl7t . 

For the range of temperatures T < M s and chemical po- 
tentials H M s that we are interested here, we find 



e s (-oo) = - 



H 



(1/2/32) 



T(l- 2/3 2 



(20) 



The interactions renormalize e s (— oo) up or down changing it 
from e s (— oo) = —H (the bare energy) of the first order of ap- 
proximation, the noninetarcting limit, to e s (— oo) = ~2(3 2 H 
when interactions are present ("dressed" energy). The renor- 
malization occurs due to interactions with other particles, 
present in the presence of large H. Therefore the asymptotic 
behavior of the free energy for the above case of H ^> M s > 
T will be 



f(H) 



[3 2 H 2 



TtT 2 
12 



(3 2 TH 
2tt 



11 
2tt 



2f3 z H 



(21) 

which leads to a constant susceptibility up to higher temper- 
atures. At H ^> M s > T (putting back also vp) the charge 
susceptibility approaches the constant value of \ — > 2/3 2 /ttvf 



X (H/A » 1) 



2(3 2 
ttvf 



(22) 



This corresponds to the limit of large density of solitons. 
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FIG. 4: The charge susceptibility of the SG model as function of the 
external field V = 2H, where H is the chemical field for solitons 
and antisolitons, is presented. The Luttinger liquid description seems 
to be not useful even at temperatures T = M s . The asymptotic 
behavior H/T S> 1 found numerically agrees with the analytical 
result x = "2/3 2 /-k, of Eq. (gl). Here (3 2 = 1/3. 
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FIG. 5: The charge susceptibility of the SG model as function of the 
external field V = 2H for the case of f3 2 = 1/2, is presented. The 
Luttinger liquid description in this case starts to break down at higher 
temperatures compared with the case of /3 2 = 1/3. 



2. Numerical results for the charge susceptibility at finite 
temperatures 

At values of the external field on the order of the gap or 
smaller and at finite temperatures the breathers and antisoli- 
tons give significant contributions to the thermodynamics of 
the system. For 1 — 1/2 for the free energy in the region 
H < M s we find 

2T(3/2) ^ (-1)" (nH\ (nM s \ 

n=l N / N / 

(23) 

where K\ is the modified Bessel function. 

Plots of the susceptibility in the whole range of the exter- 
nal field and for various temperatures compared to the gap 
are given in Figs. [4] and [5] for values of the coupling constant 
(3 2 = 1/3 and f3 2 = 1 /2, respectively. In both cases at T = 
a square root singularity occurs at V/ A = 1. This is the re- 
sult of the square root dependence of the Fermi velocity on 



the external field in the metallic phase at V/ A > 1. At finite 
temperatures the divergence in \ at the gap crossing is sup- 
pressed. At higher values of the field H however the lines for 
X obtained at various temperatures merge together. 

The more breathers in the spectrum, the faster the LL de- 
scription breaks-down. Since in the case of Fig.|5]in the spec- 
trum there are no breathers the LL description is expected to 
be a good approximation till higher temperatures than the case 
of smaller (3 2 . The comparison between Fig. [4] and Fig. [5] 
makes clear this intuitive expectation. 

At values of the chemical potential reasonably close to the 
mass gap an approximate value of the temperature for the va- 
lidity of the LL description of the model at finite temperatures 
for the case of /3 2 = 1/3 would be T - 0.1M S - 0.75K. In 
cases of stronger interactions this temperature is expected to 
become progressively lower. 
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B. Drude weight and charge-charge correlation function in I 
phase of the sine-Gordon model at T = 

Here we calculate the Drude weight and the charge sus- 
ceptibility for the gapless modes above the soliton condensate 
when the system can be described as a LL. At low energies 
and at T = the excitations of the system can be described 
by the LL action 

S = A J dtdzf-L (d t <p f - v F (d x <S>) 2 } , (24) 

where vp is the Fermi velocity in the incommensurate phase 
and K the LL parameter. This kind of LL, which is obtained 
by linearizing the spectrum around the two Fermi points of the 
soliton sea, has parameters that clearly should depend on the 
values of H and the coupling constant (3 2 of the underlining 
SG model. For the relationship of vp and K to 1 and H we 
use the exact solution of SG. 

The Drude weight and charge susceptibility gives an way 
to experimentally measure the combinations vpK and K/vp. 
They are both functions of 1 and H and we examine in this 
section this dependence. We show that the Fermi velocity in 
the I phase converges to vf very quickly with increasing H. 
vf is the Fermi velocity contained in Eq. the velocity prior 
to the gap opening. The latter parameter is the one that is 
needed for the evaluation of the gap in Eq. @. The other re- 
maining parameter needed to find M s is the value of the cou- 
pling constant 1 . As it was shown in Ref.llqthe LL parame- 
ter K converges fast to its asymptotic value K — ► (3 2 j (4/7r) as 
H increases above H c = M 8 . 

The optical absorption is proportional to the real part of the 
optical conductivity. The real part of the optical conductivity 
on the other hand is related to the imaginary part of the time 
ordered current-current correlation function by 



$te{a(uj, q x )} 



TL(iu) — > to + i0 + , q x ) 



(25) 



where the Fourier transform of the time ordered current- 
current correlation function is given by 



IL(iuj,q x ) 



dxdre 1 ' 



T T (j(x,T)j(0,0)) . (26) 



Since the electrical current can be expressed in terms of the 
right and left chiral currents by j(x. t) = vp[J(z) — J{z)] = 
(i/y/ir)d T $>(x, t), Eq. (EB! takes the form 

Re{a(q x ,uj)} = 3?m( / drdx (27) 

I nu J 

T T (d T $(x, r)<9 T $(0, 0)) e^e-^l . 

The units of the conductivity (e 2 /h) are not written out ex- 
plicitly. Eq. (|27T i leads to 



v K r 
^e{a(q x , lu)} = — Sm 

7T I 



uj + i0+ 



(lu + z0+)2 - (vp) 2 q x 



i,) 2 0? J 



(28) 
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FIG. 6: The Drude weight in the incommensurate (metallic) phase of 
the SG model as function of f3 2 and external field is shown. 



and 



5?c 



(29) 



The Drude weight will be D = vpK, where K is the Lut- 
tinger liquid parameter. 

The external field couples to the charge p(x, t) — J(z) + 
J(z) — —d x ^(x ) t)/^/TT. The charge susceptibility will be 
given by a similar equation to Eq. d27l ), now through the time 
ordered charge-charge correlation function 



drdx 



(30) 



T T {d x ^(x,T)d x ^0))e^ T e 



IUT e —1q x x 



and 



{uo + iQ+) 2 + {v F ) 2 ql 



K 

X(q x ,uj-> 0) = — T 



(31) 



(32) 



One can reach also at the same result through the thermody- 
namic Bethe ansatz equations at T = 0. The susceptibility 
can be obtained as a second derivative of the ground state en- 
ergy Eq. (l65T l with the external field. For related work see also 
Refs. 15.23.24l! 



1. Fermi velocity and Luttinger liquid parameter in gapped phase 
in the presence of chemical potential 

In the following we give the behavior of the Drude weight 
and charge susceptibility Eq. (f32t in terms of /3 2 and H ex- 
ploiting the Bethe ansatz equations for vp and K. The Fermi 
velocity for solitons in the I phase (see Fig. |2]i can be found 
by taking the derivative of the soliton branch of the spectrum 
e s (P) with P, the momentum, 



de s (P) 



dP 



de s (6) 



1 



88 dP/dO 



(33) 
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which is defined 2 ^ as 



FIG. 7: On the left panel the ratio of the Fermi velocity of I (metallic) 
phase as compared to that of the C (insulating) phase prior to the gap 
opening, v f /vf = de€s{6)\o=B /2-Kp s (B) of Eq. ([34}, is shown. 
On right panel the dressed charge (LL parameter) as function of cou- 
pling constant, or equivalently, the interedge interaction strength /3 2 
and external field is shown. The LL parameter is K — 1 resembling 
free fermions at the phase transition point when H — M s and takes 
the asymptotic value of K = 2/3 2 at larger values of H. 



de s (9) 



06 



-_ b 2tt Ps (B) 



(34) 



where Vp is the Fermi velocity of the commensurate phase 
and p s (B) is the soliton density at the Fermi point. 

The equation for the energy spectrum is given in previous 
section by Eq. dT2b and therefore de s /dO is given by 

e' s (8) + f d8'G(6 - = M s sinh(<9) , (35) 

J-B 

whereas the particle density equation is the following 

r-B 



Ps(0) 



&9'G{9 - d')p s {9') = M a cosh(6»)/27r . (36) 



<Mg) 

dH 



and as a result it is given by the integral equation 



C(0)+ / d9'G(9-6')((9') = l , 

-B 



(38) 



(39) 



— B 



K = C 2 (B)- This parameter characterizes the strength of the 
interactions of the excitations, and gives complete informa- 
tion on the exponents of various correlation functions in the 
metallic phase. The dependence of the LL parameter on the 
value of the external field is enforced by the field dependence 
of the Fermi boundary point B(H). At the insulator-metallic 
phase transition point when H = AI S the integration interval 
shrinks to a point and K = ( 2 (B — > 0) = 1. The physical 
origin of this behavior is that at this point the density of soli- 
tonic particles of the theory is very low and they behave as 
free ones, leading thus to the characteristic K = 1 value of 
the LL parameter. 

At H/M s — > +oo the Wiener-Hopf technique shows that 
the LL parameter K — ► 2f3 2 . This coincides with v F — > vf- 
Of course the limit of large values of the external field corre- 
sponds to the case of large density of solitons, and the system 
restores to properties prior to its gap opening. The cos opera- 
tor does not affect the properties of the system at these values 
of the external field. 



2. Asymptotic behavior ofDrude weight and charge-charge 
correlations 



The Fermi points ±B are defined by e(±B) = 0. We have 
solved these integral equations numerically and the graph for 
the Fermi velocity in the I (metallic) phase as a function of 
the coupling constant and external field is given in Fig.|7](left) 
whereas the Luttinger liquid parameter K in the whole range 
of H is shown in Fig. |7](right). 

From the above integral equations it is easy to obtain the 
asymptotic behaviors H — > H c and H/H c ^> 1, where H c = 
M 8 . The behavior of the soliton Fermi velocity in the I phase 
is given by 



1/2 



1 - 



4v^G(0) 



H 



- i 



(2G 2 (0) + 1) ( f- - 1 



H 



1/2 



(37) 



Notice here that the square root behavior of the Fermi ve- 
locity vp is a generic feature of all models with a gap in the 
spectrum whose low energy behavior fulfills Eq. ( fTTT l. This 
leads to other generic properties in the behavior of the Drude 
weight and charge susceptibility below. 

The parameter K is given as the square of the dressed 
charge of the solitons ((6) around the Fermi points 6 = ±B, 



I. (H/H c — > 1). Using the results of Ref.[l6| for the depen- 
dence of the LL parameter K on H at T = for H — > H c , 
the Drude weight behaves as 



D 



v F \ H, 



H 



1/2 



4V2G(0) / H 
3 \H r 



1/2 



(2G 2 (0) + 1 + 4V2 



H 



■ (40) 



The Drude weight close to the critical point H c increases as a 
square root function of (H/H c — 1). The same holds for the 
Fermi velocity in the I phase, v F , close to H c . 

At T = 0, the H — > H c asymptotic behavior for the charge 
susceptibility is 



y/2: 



-1/2 



4%/2G(0) / H 



(41) 



1/2 



1 - 4V2G(0) 



50G 2 (0) 
9 



£-1 
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FIG. 8: The soliton-antisoliton contribution to the optical conductiv- 
ity in the C phase for the SG model. 




FIG. 9: The breather-breather (e\€2) contribution to the optical con- 
ductivity in the C phase for the SG model. 



which clearly shows that it diverges as a square root function 
of (H/H c —1). At small values of the coupling constant 1 we 
expect the square root behavior for D and \ t0 hold only for 
a very small interval of H above H c . For larger values of H 
this behavior is expected to change to a logarithmic function 
of (H/H D - 1). 

II. (H/ M s > 1). At H /M s > 1 it was shown previously 
in Ref. 16 that the LL parameter K — > 2(3 2 . Therefore at 
H ^S> M s we expect the experiments on the Drude weight and 
the charge susceptibility to give the following combinations of 
the coupling constant and Fermi velocity prior to gap opening 

2d 2 

D = 2{3 2 v F , X = — , for H > M s . (42) 
ttvf 

(Here the usual 87r of Ref.llqis absorbed in (3 2 , which belongs 
now to < (3 2 < 1). Measuring therefore the Drude weight 
and the charge susceptibility at large values of the external 
potential should define the quantities needed for measuring 
the mass gap of the experiments of Kang et alX. 



C. Optical conductivity in the insulating phase 

One important quantity to study experimentally in the con- 
text of the applicability of the sine-Gordon model in QHLJ 
systems is the optical conductivity. This quantity is quite use- 



ful since one can extract valuable information about the polar- 
ized excitations of the system. 

Just as in the metallic phase of the previous section, the 
optical conductivity in the insulating phase is defined as the 
Fourier transform of the current-current correlation function. 
Here one has to use the integrability property of the SG model, 
the form factor approach 32 . The particles that contribute will 
be the ones that are odd under charge conjugation, the odd 
breathers e^ n —\ in the SG model. The single breather contri- 
bution is simply given by a delta-peak at a frequency given by 
the mass of the breather (Mk = 2M s sm(^£k)). Other contri- 
butions will be the two-particle ones like, soliton-antisoliton 
and the odd-even breathers. At higher energies there will be 
contributions also form higher multiparticle states, like the 4- 
particle ones etc. 

The optical conductivity for the sine-Gordon model is cal- 
culated by Controzzi, Essler and Tsvelik 31 . We use here their 
results to illustrate plots of the first-second breathers and the 
soliton-antisoliton contributions to the optical conductivity as 
shown in Figs. [8] and [9] More detailed calculations for higher 
order terms are given in the appendix. 

Let us take as an example the case of f3 2 = 1/4, where 
there are two breathers in the spectrum. As function of the 
frequency, the zero temperature optical conductivity will be 
zero for u < M\ = 2M S sin(£7r/2), a delta-function peak at 
u) = Mi and then a continuous weight at uj > 2M S for the 
soliton-antisoliton which comes before the contribution from 
the first-second breather u > M\ + M 2 = (1 + V3)M S . 

At lower values of j3 2 (stronger interedge forward in- 
teractions in the junction), the mass of the breathers de- 
crease whereas the mass of the solitons increases, the num- 
ber of breathers on the other hand increases and the odd-even 
breathers contributions will arise as continuum contributions 
until the soliton-antisoliton continuum contribution. Due to 
the heavier mass, the latter contribution in the optical con- 
ductivity will be much stronger than the one of the odd-even 
breathers. In the experiments, optical absorption at small fre- 
quencies will signify the presence of breathers in the spectrum 
and thus signify the presence of strong interedge Coulomb in- 
teractions in the line junction. 



D. QHLJ at filling fraction v ~ 2 

The QHLJ at filling fraction v ~ 2 in the case of the pres- 
ence of a single impurity along the barrier has been discussed 
beautifully by Kim and Fradkin in Ref.lSJ Here we are mainly 
interested in the case when interedge tunneling can take place 
at any point along the barrier but largely follow the ideas of 
Ref.ll 

The calculations of the previous sections can be extended 
also to the spin-full case, when the 2DEG subsystems on the 
sides of the barrier are at filling fractions v ~ 2. The Hamilto- 
nian in these ID systems separates into two commuting parts, 
the spin and charge ones H = H c + Ti s - We focus here on 
the case of a gapful spin sector, where the chemical potential 
for the topological charges consists naturally on the Zeeman 
coupling term. Information on the strength of electron inter- 
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actions can also be obtained by working in this sector. We 
discuss here the possibility for observation of the SILL phase 
that has been found to exist in ID systems with spin-charge 
energy scale separation^. 

The edge states of a 2DEG on a QH plateau at bulk filling 
v ~ 2 are composed of spin up and spin down branches. To 
account for both of them we need to introduce two bosonic 
fields, (j) a , where a =1,1 are the up and down spin pro- 
jections. The right and the left movers are represented by 
R a ~ e-'v 7 ^-, and L a ~ e Jv/3¥< K respectively. The elec- 
tron spin brings about changes in the dynamics — terms like 
the spin-spin exchange interactions as well as the Zeeman in- 
teractions have to be added into the Hamiltonian. Introduc- 
ing the scalar U(l) and SU(2) vector currents J =: R) a R a :, 
J =: R) a (ij a p/2)Rp :, and similarly for the left movers, 
R — ► L, the charge and spin components of the Hamiltonian 
(without the Zeeman term) can be written as 



n s = 



2 
2ttv 



■(J J + J J + g c JJ) 



(33 + JJ + g s 33) 



(43) 



(44) 



In the standard way (Ref. 27) we introduce the charge and 
the spin bosonic fields tpi = (<pf ± <pi)/\/2 for i = c,s and 
similarly for the opposite chirality fields fa = (fa ± <j>y ) / 
The charge excitations remain free and described by the LL 
model. 



H r = 



dx 



K c {d x <S> c ) 2 + -^(d x Q c ) 2 



(45) 



where $ c = <fi c + (j> c and 6 C = cf) c — fa. The Fermi velocity 
is renormalized by the g c density-density type interactions. 
(The relationship between v c , K c and g c is nonuniversal — 
when the interactions are weak these parameters are related 
by simple expressions (Ref. l27l)l . 

In H s the g s type perturbation is subject to renormalization 
group (RG) flow. JJ is a marginal operator and the sign of 
interactions is important in that it can flow the system toward 
the strong coupling phase, with no scale invariance, or it can 
scale it to zero (the gapless phase). In general the spin-spin 
exchange interaction g s 33 can be extended to an anisotropic 
one 



,g s JJ -> g\\J z J z 



1 



;9± (J + J~ 



h. 



This can be the case in the samples of the experiments of Rang 
et alX. The new constants are: g± = g s the XY component 
of exchange interaction, and g\\ = (1 + X)g s the Ising type 
coupling. A is the strength of the Ising anisotropy. The mag- 
netic field introduces a Zeeman type coupling [Eq. [9). The 
spin part of the Hamiltonian is described by the sine-Gordon 
model, where the cos operator is introduced by the XY-type 
spin-flip backscattering interaction g±. The RG flow of the 
model is discussed in detail in Ref. 9. The isotropic case 
(A = 0) is subject to SU(2) symmetric RG flow, with pertur- 
bation marginally irrelevant at the ferromagnetic case g s < 0, 
and marginally relevant in the antiferromagnetic case g s > 0. 



The latter type of interaction generates a gap in the spin sec- 
tor. A gap in the spin sector can open up even in the case of 
ferromagnetic type interactions, when the magnetic interac- 
tions are anisotropic. This happens when A < 0. Either of 
these cases would be of interest in the context of our previous 
discussion. 

The Zeeman coupling (magnetic field) can be regarded as 
a field that is coupled with the topological charges of the SG 
model 



Ax 



K s {d x <$> s ) 2 



1 



(d x e s y< 



(47) 



-2/x s cos(v87r$ s ) , 



We discuss here various cases that arise concerning electron 
tunneling conductance across the junction. 

In systems where the electron interactions are very strong, 
a separation in the energy scales in spin and charge sec- 
tors arises. At low temperatures it is possible for the en- 
ergy to be low compared with the characteristic charge en- 
ergy and high compared with the characteristic spin energy, 
E s <C ksT <C E c . In this window of energy, the charge sec- 
tor is very close to the ground state whereas the spin sector is 
thermally excited. As it was shown in Ref. [lj this state has 
quite interesting properties, particularly the correlation func- 
tions decay exponentially in space and the exponents of the 
tunneling density of states is also different from a simple Lut- 
tinger liquid. In the QHLJ systems that we examine here, 
such conditions can arise naturally. Due to the small barrier 
widths, less than a magnetic length Ib in the systems of Ref.Ql 
we believe that Coulomb interactions are strong whereas the 
energy scale in the spin sector can be varied by tuning the 
magnetic field. The charge energy scale is of the order of 
E c ~ e 2 /els ~ 10.76meV. On the other hand, the spin 
excitations in the case of the antiferromagnetic g s > inter- 
actions are gapped. At values of the Zeeman term bigger than 
a critical value H > H s , the gap closes. It is at values of 
the Zeeman field only slightly in excess of Q c that the SILL 
phase can occur. At H — > H s from above, the characteris- 
tic Fermi velocity of excitations in spin sector becomes very 
small, v s ~ vp[(H/H s — l) 1 / 2 + . . .] ^0 and the energy 
scale in spin sector vanishes E s ~ %v s /Ib — * 0. 

As discussed in Ref. 9, the most relevant tunneling opera- 
tors that account for both charge and spin transport across the 



(46) junction are the electron tunneling O c \ = R1l 



r\l 



h.c. ~ cos(\/27r K c <&c) cos(v / 27rKj$ s ), the singlet pair tun- 
neling dp a i r = i?ji?jL|L T + h.c. — cos(V87ri^ c $ c ), 
and the triplet pair tunneling O s = i?j if ij-Rj. + h.c. ~ 
cos( v / 87riir s $ s ). 

We analyze in the following the relevance of the tunneling 
operators at various values of the Zeeman term in the spin 
gapped state. Since the spatial separation of the counterprop- 
agating edges on the sides of the barrier is of the order of 
one magnetic length or less we assume interedge tunneling 
can take place at any point along the barrier. The scaling di- 
mensions of the operators depend on the LL parameters of the 
charge and spin sectors and are given by 9 d c \ = (K c + K s )/2, 



12 



whereas d pa ir = 2K C and d s = 2K S . First of all, in the 
spin gapped phase the only allowed tunneling process is the 
singlet pair tunneling O pa i r which is independent of the spin 
bosonic field. The spin gap is generated by the operator 
cos(v / 87i r K7<i> s ) [either when g s < and any A, or when 
g s > and A < 0], and therefore <1> S takes as value one of the 
minima of the potential. The average of the cos(v / 2lrK7$ s ) 
will alternate in this set and result in a vanishing average. 
Opair term will lead to a gap opening in the charge sector 
whenever K c < 1. The charge and spin currents will be trans- 
mitted across the barrier. By sweeping the Zeeman term to 
higher values than the spin gap in the region H ~ H s , at 
finite temperatures, due to vanishing energy scale in spin sec- 
tor, the SILL should arise. Due to thermal excitation the spin 
correlation functions here again decay exponentially and spin 
tunneling processes are suppressed. The only allowed tunnel- 
ing processes will be again the singlet pair tunneling although 
this will be relevant only at K c < 1/2. The jump in scaling 
dimensions of this gap generating operator in charge sector is 
a signature of the SILL state. A perfect charge transport at 
these values of the Zeeman term will signify strong correla- 
tions along the barrier. At higher values of the magnetic field 
the pair tunneling operator and the electron tunneling operator 
will be both relevant. At strong interactions we expect O pa ; r 
to be more relevant then O c \. Charge and spin transport will 
display a peak as function of the bias voltage. 

In the presence of a single impurity the phase diagram as 
function of K c , K s , and magnetic anisotropy, was studied in 
Ref.Hl Here we add that in the spin gapped phase in the pres- 
ence of a single impurity and the SILL phase (H ~ H s ) a 
divergence on the density of states at the impurity point arises 
whenever K c > 1/4 and charge tunneling through the barrier 
at low bias voltages does not occur. In the opposite case at 
K c < 1/4 the presence of the impurity will lead to a suppres- 
sion of the density of states and a peak in the charge conduc- 
tance through the barrier should be observed. 

IV. CONCLUSIONS 

In this paper we studied the properties of QHLJ created by 
the narrow cleaved edge overgrowth barriers. Such line junc- 
tion have a sort of properties that make them be advantageous 
to other sorts of line junctions, as for instance, the gate created 
junctions, in that, first, the barriers present to the electrons suf- 
ficiently abrupt potentials that edge reconstruction phenomena 
do not occur. This allows for the edge states of the CEO cre- 
ated line junctions to be described by single chiral channels. 
Second, these barriers can be narrow enough to allow for elec- 
tron tunneling to take place at any point along the length of the 
barrier, leading to low energy description in terms of general- 
ized sine-Gordon models. And third, these systems allow for 
strong Coulomb interactions that lead to pronounced features 
in transport experiments. 

We discussed here ways of estimating the parameters of the 
field theory model that describes the QHLJ system 1 at low en- 
ergies. These parameters are hard to estimate with accuracy 
from microscopic models. Therefore we undertook an indirect 



way for their estimation. We showed that information on the 
value of the LL parameter and the Fermi velocity can be ex- 
tracted by using the integrability property of the sine-Gordon 
model. When forward scattering repulsions are present the SG 
model is characterized by a gap in the spectrum. The excita- 
tions consist on topological particles and their bound states (if 
Coulomb interactions are strong). The gap can close in the 
presence of fields that couple with these charges. In the so 
obtained metallic phase we propose the study of the charge 
susceptibility and the Drude weight. It turns out (see Figs. [6] 
Q that not very far from the critical value of the gap clos- 
ing, the LL parameter which depends on H and 0, reaches as 
its asymptotic value, the one that would characterize the sys- 
tem prior to gap opening. In the insulating phase, to obtain 
information on the properties of the system, we propose the 
measurement of the optical conductivity of the line junction. 
The optical conductivity provides valuable information on the 
polarized excitation of the system. Delta-function type fea- 
tures appear in optical conductivity for each particle type odd 
under charge conjugation (all odd breathers e2 n -i), as well 
as two and more particle bands appear in the spectrum. This 
measurement should therefore lead to information about the 
number of particles present, and thus the strength of the inter- 
actions. Here the form factor formalism 32 is used. 

We discussed here also the case when spin degrees of free- 
dom take part in the dynamics, i.e. the case of bulk-filling 
fraction v ~ 2 on the 2DEG subsystems. Most interesting is 
the study of the case a spin gapped phase. In this case the Zee- 
man term can be tuned to values just above the mass gap, re- 
sulting in a strong separation of charge and spin energy scales. 
At finite temperatures, much less than the charge energy scale 
but larger than the spin energy scale, the SILL state can be 
realized. Such a state has been proposed to exist 12 in ID sys- 
tems with strong Coulomb interactions - Wigner crystals in 
charge sector and weak spin exchange interactions, but its ob- 
servation is yet to be seen. Its unusual properties stem from 
the reduction of the spin energy scale, which turns out to have 
important effects on the correlation functions and thus trans- 
port properties of this state. 

We focused here on the coherent line junctions where tun- 
neling takes place at all points along the junction. At wider 
barriers, due to momentum conservation, interedge tunneling 
is prohibited except at points of impurities along the barrier. 
The low energy physics of systems with just a single impurity 
can be mapped to the boundary sine-Gordon model (Ref. |9[), 
which allows in principle for an exact calculation of the dif- 
ferential conductance. In real experiments however, a tem- 
perature and voltage dependence of the tunneling amplitude 
occurs. The tunneling conductance calculated with the Bethe 
ansatz taking into account these other effects will be discussed 
in a future publication 3 ^. 



ACKNOWLEDGEMENTS 

The author is grateful to Dave Allen, Matthew Grayson, 
Woowon Kang, Allan MacDonald, Stefano Roddaro, Tilo 
Stroh, and Alexei Tsvelik for illuminating conversations. The 



13 



author is especially grateful to Allan MacDonald for teaching 
him the Hartree-Fock theory in the context of line junctions 
in QH systems and many other illuminating discussions and 
a careful reading of the manuscript and to Eun-Ah Kim for 
valuable comments on the manuscript and discussions. This 
work was supported by the National Science Foundation un- 
der Grant No. DMR-0412956. 



V. APPENDIX 



TBA equations for purely elastic theories and the T = 
limit for the sine-Gordon model 



Here we write the thermodynamic Bethe-Ansatz (TBA) 
equations for systems for which there is no reflection on scat- 
tering of the particles of their spectrum, i.e. we consider only 
diagonal scattering or purely elastic theories. These equations 
apply for the SG model only at particular values of the cou- 
pling constant, namely (3 2 = 1/n, where n is an integer or, in 
the presence of the external fields, only at temperatures much 
lower than the size of the gap, T <C M s . At zero temperature, 
which is the case of our interest in (uTO . the BA equations 
for SG in the presence of external fields can be derived only 
based on the fact that in case H > h c = M s in the ground 
state there will appear a soliton condensate and therefore the 
equations with a single kind of particles in the spectrum ap- 
ply. The T — will be good for every value of the coupling 
constant of SG. 

We are dealing here with a Lorentz invariant model and in 
this case it is usual to use the following parametrization for 
the particles 's energy and momenta 



E = m cosh 9 



P = (tti/vf) smh9 



(48) 



where m is the mass of the particle and vp the Fermi ve- 
locity [corresponding to the Fermi velocity of the C phase 
of Eq. (@)]. In this parametrization Lorentz invariant quan- 
tities like the scattering matrix elements will depend only on 
the differences of rapidities (invariant under Lorentz boosts 
which shift the rapidities by the same constant). Suppose we 
have N particles of mass m on a closed circle of length L. Let 
the two-body scattering matrix of particles with momenta pi 
and p2 be equal to S(9\ — 82) — exp{i</?(6>i — #2)}- If we 
consider a particle of momentum pj traveling around the cir- 
cle, on returning to the initial point its wave function acquires 
the factor 



exp 



(49) 



and this has to be equal to unity. Therefore the argument of 
the exponential has to fulfill 



rriL sinh^i + <j>(0i — 6j) = 2-Ktii 



(50) 



(51) 



There are N different integer numbers Hi (eigennumbers) 
which determine the eigenvalues of the Hamiltonian 



£ = ^3 m cosh t 



(52) 



In the thermodynamic limit L — ► 00, the distance between the 
rapidities, solutions of Eq. (15 It . is of the order of 1/L and in 
this case one introduces the densities of rapidities occupied 
by particles p p (8) and the density of holes ph{8) (rapidities 
which can potentially be occupied by particles). Taking the 
derivative of ( BTT l with 8i one gets the following integral equa- 
tions for the dressed distributions 

Pp (9)+p h (6) = 1 ^^- + (K*p p )(d) , (53) 

Z7T 

and for the more general case of several types of particles a 
with a = 1,2, ... ,n, with masses m a 



6=1 



where (m a cosh 9) /2vr = (l/2vr)dp°/d6> = p° a:P (9) is the 
bare distribution function of the particles of type 'a', whereas 
the Kernel K is related to the two-body scattering matrix 



(55) 



The energy of the system is given by 

£ = I ^2 p%(0)m a cosh 6 d6 , (56) 

"* a 

and in the presence of a chemical potential coupled to the 
"charge" of the particles it is 



£ = £ - E ha<la 



(57) 



In Eq. d56i l m a cosh (9 = cq(9), the bare energy. 

The entropy of the distribution between 9 and 9 + d9 is the 
logarithm of the number of distributions 



dS = In 



{L[p a p (9) + p a h (9)}d9}\ 
[L Pi {9)]\[Lpl(9)]\ 



(58) 



Here one uses the Stirling formula ln(n!) ~ n(ln n~ 1). Then 
the entropy per unit length is 

S / L = I £ { W + ^ [ P a p (0) + p a h (0)] 

J a 

-p a p {9)\np a p {9)-pl{9)\np a h {9)}A9 . (59) 

From the condition of the minimum of the free energy F = 
E—TS in the equilibrium state one gets relations between Sp p 
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and Spfr, which when substituted in ( 154-b lead to the following 
set of integral equations 



m a cosh9-h a q a = e a (9)+T^2 K ab *\n (l 



6=1 



+ e 



-e b /T 



where the functions e a (0) are defined by 
= cxp j — 



(60) 



Using the dressed energies one finds a formula for the free 
energy per unit length 



F/L = -Tj2m a J — cosh 6 In 



d6 



l + e 



(62) 



As T — > the equations for the dressed energies (f60l > reduce 



to 



i a cosh 9- h a q a = e a {9) - if afc * e b , (63) 



6=1 



where e~(6) is the negative part of the dressed energy e a (0) 
[usually between the interval we note (— B, B)]. If from the 
set of the particles only one has positive charge, then all the 
energies except ei(0) are positive everywhere and the set of 
the equations d63l reduces further to 

,B 

ei(6)- d9'K 11 (9-9')e 1 (9') = m 1 cosh9-h 1 q 1 , 

J-B 

(64) 

where ei (±-B) = 0. From Eq. d62l one finds the soliton (par- 
ticle a = 1) contribution to the ground state energy of the 
system 



£ = / A9p° lp {8)e 1 {9)=m 

—B 



66 



i / — cosh6>ei(6>) 

—B 27T 



(65) 

In Eq. d65b . p® p (9) is the bare distribution of particles of type 
'1' in the ground state. 

In the following we give the 2-particle soliton-soliton and 
soliton-antisoliton scattering matrices for the SG mode l 18 ! 19 



S s a s s (6) = Sll(6) = S (8) , 

sinh(0/£) 



sim = sim 



sinh | ( 



sinh(wr/£) 
sinh 4(0 — iir) 



(66) 

-5 O (0) Trans. 
I 

S (9) Refl. 



with 



S o (0) =-exp 



sin(0i/7r) sinh(i^t) 
tcosh(i/2)sinh(£t/2) 



At £ = 1/n the Sjj(0) = Sjf (0) = and the scattering is 
purely transmissive. In the above equations 



il 2 



1-p 2 



(67) 



At 1 = 1/2, which corresponds to the fermion free point, 

S ss = -1 and S sS = S = -1. 



(61) B. Kernels of particle interactions at 1 — 1/3 and /3 2 = 1/5 



To calculate the interacting kernels of the particles we use 
the following formulas (Ref.[33l) 



sinh(n£) 



[(n-l)/2] 



fc=0 



sinh£ 
and (Ref.|33| 

d£ cos a£ 



(-l) fe ("-^ 1 )2"- 2fc - 1 cosh"- 2 ^ 1 £ , 



cosh/3^ _ 7T cos cosh 2£ 



cosh 7^ 7 C0S hS£+cos 



CiTT 



(68) 



, (69) 



for [ | Re/3 1 < Re7, for real a]. 

For the 1 = 1 /3, the interparticle interacting kernels are 



K ss {9) 



K bs {9) = 
1 1 



2tt cosh 9 



\f2 cosh6» 
Vcosh(20) ' 

K bb {6) = - ' 1 



7r cosh 



(70) 
(71) 



For the 1 = 1/5, the interparticle interacting kernels are 



K ss (9) 
K ls (0) 
K 2s (9) 
K 3s (6) 
K n (9) 

K 22 {9) 
K 33 (9) 



1 / 1 



2V2cosh<9\ 
'tt ^ cosh 9 + cosh(20) J 

1 / 4cos(7r/8)cosh6» \ 
tt Vcos(tt/4) +cosh(2<9) / 

1 / 2 2\/2cosh6>A 
'tt I cosh 6» + cosh(20) J 



1 / 2 % /2cosh6»\ 
^\ cosh(26») J ' 

2 / 1 2V2cosh6» N 
~^ 1^ cosh (9 + cosh(20) 

1 / 4 6V2cosh6l' 
~^ I cosh 9 + cosh(20) 



(72) 
(73) 
(74) 
(75) 
(76) 

(77) 

(78) 
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C. Optical conductivity of SG model in the gapful regime 

In the massless regime we computed physical quantities 
like the Drude weight and the optical conductivity. Similar 
calculations can also be done in the massive regime. Here the 
integrability property of the sine-Gordon model can be used. 
Indeed, in addition to the energy spectrum that is obtained 
through the Bethe ansatz, quantities like the form factors and 
correlation functions can be computed exactly based on the 
bootstrap formalism. In this section, the low frequency opti- 
cal conductivity is computed exactly using the form factors. 
Similar calculation has been done previously 31 to determine 
the optical conductivity of the Mott-Hubbard insulator. We 
try in this appendix to extend their work to compute additional 
breather-breather contribution. 

The optical conductivity is defined as the Fourier transform 
of the current-current correlation function. It is an important 
physical quantity that provides information on the polarized 
excitations of a system. Explicitly, the relationship between 
the optical conductivity and the current-current correlation 
function is given by: 



1 ( f + °° 
— Im( i I dx 

U) 



x\j{x,t), 3(0,0)} 



(79) 



The current-current correlation function, can be computed by 
introducing the identity resolution, as follows 



n— a 



l(0|j(0,0)|( 



1 ■ ■ ■ u n/ai ...an 



(80) 



where |0i . . ■0 n )ai...a„ is a n particles state with rapidities 
6*i . . . 0„. The indices a, take for value ±1/2 corresponding 
to a soliton or anti-soliton or a, = 1, . . . , [l/£] corresponding 
to the different breathers, pj and ej are, respectively the mo- 
mentum and energy of the jth particle. The link between these 
variables and the rapidity is given by Eq. (42). The matrix ele- 
ment (0|j(0, 0)|#i . . . n )oi...an i s called form factor and will 
be noted f ai ... an (0i • • ■ 9 n ). The knowledge of all form fac- 
tors leads to the current-current correlation function and thus 
to the optical conductivity. From the previous definitions, one 
easily deduces the following relation: 



\ n=0 at J y 



A6 n 



fai- ■ ■ a n 



lo + ie — J2j Mj cosh 8j 

S(q + J2j Mj sinhOj/vp) 
to + ie + V ■ Mj cosh 9 j 



(81) 



In the limit e — > 0, the imaginary part can be computed lead- 
ing to an additional Dirac function over the frequency. The 
one-particle and two-particle states contribution to the con- 
ductivity at vanishing momentum is given by: 



c 



TTV 

M 2 



F |/ a (0 = O)| 2 <5(u;-M a ) (82) 



v- 2ttv$ \f aua2 (6 1 ,e 2 )\ 2 e(u; - Mi - M 2 ) 
^ lo [(a; 2 - (Mi - M 2 ) 2 ) (w 2 - (Mi + M 2 ) 2 )] 



The particle rapidities in the two particle contributions are re- 
lated to each other and to the frequency by the following rela- 
tions: 



Misinhfli + M 2 sinh6> 2 = 
Micosh#i + M2Cosh#2 = ^ 



(83) 
(84) 



The limitation of the computation to the one particle and two- 
particles states is valid at least at low frequency, u. More 
explicitly, here are some rules to determine the form factors 
required for the computation of the conductivity: 

• Only the form factors, (0|j(0, 0)|0i . . . 9 n ) ai ... an , 
where the sum of the n particles mass is less than lo 
will contribute to a(u>, q). 

• The current, j(x,t), being odd under charge conjuga- 
tion while the ground state, |0) is even, only the states, 
(0|j(0, 0)| 0i . . . 0„)ai-o n , odd under charge conjuga- 
tion will lead to a non-vanishing form factor. 

The one and two-particle states odd under charge conjugation 
are the following: 

• |01, 02)l/2,-l/2 — |01) 02)-l/2,l/2 

• \0)n where n is odd. 



where r + s is odd. 



It is not surprising that the lowest breather (n = 1) is odd 
under charge conjugation. The soliton and anti-soliton be- 
ing fermions, the multi-solitons states have to be odd under 
particle exchange. This holds also for bound states like the 
breather. It implies that if the state is odd under charge con- 
jugation, it has to be even under parity transformation and 
vice versa. However, since there is an attractive short range 
interaction between the particles, the lowest bound state (or 
breather) is obtained when the particles position overlaps the 
most. This overlap is maximized for the states that are even 
under parity and thus odd under charge conjugation. The 
soliton-antisoliton form factor has been known for a long time. 
Explicitly, it is given by: 



'|i(o,o)|0i 



M. 



72 )ss 



cosh-j- 



71" cosh- 



-F(9-) 



(85) 
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where M s is the soli ton mass, 9± = 9% ± 62, and F is called 
the minimal form factor which is given by: 



F{9) = sinh 



xexp 



(86) 



sinh 2 (t(l - i0/7r))sinh((l - £)*) dt 



sinri(£t)sinh(2i)cosh(t) 



t 



One implication of the bootstrap formalism is that all breather 
form factors can be deduced from the soliton-antisoliton form 
factors. In particular, the form factor for the kth breather (k 



J 



being odd) is given by: 



(O|j(O,O)|0) fc = 2M fl £cosh(0) 



nf=i tan 1 ! 5 



1/2 



F(ik£-iir). 



(87) 

Similarly, the two breathers form factor can be obtained from 
the four soliton-antisoliton one. The procedure for its calcula- 
tion has been described by Smirnov. Following the procedure, 
we obtained the following result (we assume that k% is odd and 
&2 even): 



<0|j(0, 0)\9 U 9 2 ) 

fel ,/C2 



(88) 



x E 

3=0 L 



nfio sinh 5 ( e - + ¥( fc i - fc 2 + 2j - 20) nr= ,^ sin w - n 



where 



fccos' 



dk = 



and 



Cfci,fc 2 (^) = Ck u k- 



x exp 



Cfci.fea = ex P { 2 



2£ 3 sin| 



1/2 






f r 


exp< 









sinh 2 ((2 - 

iinh((l-gfc)2t) 
2 



sinh(2i) 



sinh((l - £)t) dt 
cosh (i) sinh (£t) i 



2njiv fc2 ' sinh 2- (* + - fc 2i - 2i)) 



sinh 2 ((l - i0/7r)i)cosh((l - £fci)f)sinh(£M)cosh(£t) 



(1 sinh(£i)cosh(i)sinh(2t) 
fc2Cosh(i)sinh(£i) — cosh((l — £fci)i)sinh(££;2i)cosh(£i) dt 



sinh(2t)sinh(£t)cosh(i) 
I 



t 



(89) 



(90) 



(91) 



Using the obtained form factors, the different contributions the soliton-antisoliton contribution to the optical conductivity, 
to the conductivity are shown in Figs [8] and [9] Fig. [8] shows This contribution dominates at frequencies close to 2M S . 
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For temperatures T < M s and large values of the external fields, 
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After making use again of dl7t one gets 
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The first part of Eq. ( |93t can be written as 
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which is proportional to the Rogers dilogarifhm. After substitut- 
ting 
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